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Electron transport through short, phase-coherent metal-graphene-metal 
devices occurs via resonant transmission through particle-in-a-box-like states 
defined by the atomically-sharp metal leads, we study the spectrum of particle-in- 
a-box states for single- and bi-layer graphene, corresponding to massless and 
massive two-dimensional (2d) Fermions.The density of states D as a function of 
particle number n shows the expected relations D{n) ~ n V2 for massless 2d 
Fermions (electrons in single-layer graphene) and D(n) ~ constant for massive 2d 
Fermions (electrons in bi-layer graphene). The single parameters of the 
massless and massive dispersion relations are found, Fermi velocity v F = 1.1 x 
10 6 m/s and effective mass m* = 0.032 m e , where m e is the electron mass, in 
excellent agreement with theoretical expectations. 



The problem of a wavelike particle confined to a hard-walled box is one of 
the most basic problems in quantum mechanics. The spectra of the particle-in-a- 
box are strikingly different for massive and massless particles: massless particles 
(e.g. photons) have energies which depend linearly on quantum number, while 
the energies of massive particles (e.g. free electrons) depend quadratically on 
quantum number. For Fermions in two dimensions (2d) this leads to a density of 
single-particle states D = dn/dE, where E is the particle energy and n the particle 
density, which varies as the square root of particle density, i.e. D(n) ~ n V2 , for 
massless particles, and is independent of particle density for massive particles, 
i.e. D{n) ~ constant. Here we show that mesoscopic, ballistic 1 single-layer 2,3 and 
bi-layer 4 ' 5 metal-graphene-metal devices act as Fabry-Perot cavities for electrons 
confined between the atomically-sharp partially-reflective metal leads. Electronic 
conduction occurs through resonant states of the Fabry-Perot cavity, which are 
exactly analogous to the particle-in-a-box states of an electron confined by 
perfectly reflective walls. D(n) is measured, and the expected dependences on 
particle number are verified: D ~ n V2 for massless particles in single-layer 
graphene, and D ~ constant for massive particles in bi-layer graphene. D{n) is 
used to extract the single constants in the dispersion relations: the Fermi velocity 
v F = 1.09x10 6 m/s for massless particles in single-layer graphene 2,3 , and the 
effective mass m* = 0.032 m e , where m e is the electron mass, for massive 
particles in bilayer graphene 4,10 in excellent agreement with theoretical 
expectations 5,10,13,21 ' 22 and other experimental results 2 ' 3,19,20,23 . 



We first review the results of the two-dimension particle-in-a-box problem. 
Figs. 1A and 1B, illustrate the massless and massive dispersion relations 
respectively: 
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where f) = hl2u, and h is Planck's constant; each dispersion relation is 
characterized by a single parameter, m* for the massive dispersion, and v F for 
the massless dispersion. For particles confined to a 2-dimensional box of width 
W and length L the hard-wall boundary condition quantizes the wavevector 

k = (k x ,k y )=( p — ,q— resulting in two positive quantum numbers p, q. Figs. 1C 

and 1D illustrate this quantization, where each point represents an allowed 
wavevector. Then the energies in terms of quantum numbers are given by the 



familiar relations E = hv F n A 
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(massive). For Fermions at zero temperature, the occupancy of particle-in-a- 
box state will be the degeneracy of individual states g for state of energy E < E F 
(the Fermi energy) and zero for states E > E F . The number of states with E < E F 
is given by N ' = gk^WL I An , where /c F = /c(E F ) is the Fermi wavevector. Figs. 1C 
and 1D illustrate the occupied states included for equally-spaced values of E F , 
and Figs. 1E and 1F show the energies of the particle-in-a-box states as a 
function of particle number N for massless and massive 2d Fermions respectively. 
The linear and square-root dependences of E(/V) for massless and massive 2d 



Fermions respectively are evident in Figs. 1E and 1F. Using the areal density of 
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particles n = — = ^_l we then have the following relations for the dependences 

WL An 

of the Fermi energy E F and density of states D on density: 
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E F =v^ F yfm (massless) (2a) E F = (massive) , (2b) 

2m * 

D = —^— j— (massless) (3a) D = (massive). 

(3b) 

Thus the measurement of D as a function of n distinguishes massive and 
massless particles, and (given knowledge of the degeneracy g) also determines 
the constants of the dispersion relations v F and m*. 

Single- and bi-layer graphene may be used to realize the dispersion 
relations in Equations 3a and 3b as follows. Single-layer graphene is well- 
described by a tight-binding model considering only the Tr-orbitals at each atomic 
site. At zero doping, the tt and tt* bands meet at two points in the Brillouin zone 
with wavevector K. This crossing is preserved as long as the two atoms A and B 
in the unit cell are equivalent. Taking E(K) = 0, and measuring k away from the 
K point, the band structure is well-approximated by Eqn. 1a, with 

Vp =(V3/2]b7 /r| ~ 1.0 x 10 6 m/s where a = 2.46 A is the graphene lattice 

constant and y ~ 3.16 eV 9 is the nearest-neighbor hopping parameter. In 
Bernal-stacked bi-layer graphene 4 ' 5 ' 10 ' 21 ' 22 atom A in one layer is stacked above 
atom B' in the 2 nd layer, and this A-B' coupling breaks the AB equivalency of the 
graphene unit cell and results in two bands which may be approximated as 



10 , where Yi ~ 0.4 eV 24 is the inter-layer 

(A-B') hopping parameter. At k = the effective mass is given by m =yJ2v 2 F ~ 
0.03 m e . In both single- and bilayer graphene the degeneracy g = 4, due to the 
two-fold spin degeneracy and the two-fold valley degeneracy (presence of two K 
points). 

We now discuss the graphene samples used in this study. We 
mechanically exfoliate Kish graphite on 300nm Si0 2 /Si substrates to obtain 
single and bilayer graphene 2,3,4 . Single layer graphene is more transparent than 
two or more layer graphene under optical microscope as seen in Fig. 2. After 
locating graphene flakes, Cr/Au(5nm/50nm) were thermally deposited for 
electrical contacts. The channel lengths L for Fabry-Perot interference 
measurement are 200nm - 300nm and measured by scanning electron 
microscope. The maximum field-effect mobilities at low temperature estimated 
from the four -probe resistivity of the adjacent graphene sections are 15,000 
cm 2 /Vs and 4,000 cm 2 A/s for single- and bilayer graphene respectively. Figs. 2A 
and 2B show completed single- and bilayer graphene devices respectively. 
Electrodes are patterned on each graphene flake to form a large-area Hall-bar 
arrangement for characterizing the longitudinal and Hall conductivities (o xx and 
o xy ) of the sample. In addition, pairs of closely-spaced (150-300 nm) electrodes 
act as Fabry-Perot cavities on the same sample. 

Figs. 2C and 2D show o xx and o xy for the single- and bi-layer graphene 
devices shown in Figures 1a and 1b, respectively, measured in high magnetic 
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field (9 T) and as a function of back-gate voltage V g , which controls the carrier 
density n = c g V g /e, where c g = 1.1 x 10" 8 F/m, and e is the electronic charge. The 
quantized Hall effect (QHE) is evident as plateaux with o xy = ve 2 /h, and 
corresponding minima in o xx - Berry's phases of tt and 2tt lead to QHE in single- 
and bi-layer graphene at filling factors v = 4(/'+1/2) and 4(/'+1), where / is an 
integer[refs], thus our observation of the half-integer and full-integer QHE 
confirms the identification of these samples as single- and bi-layer graphene 
respectively. 

Figs. 2E and 2F show the two-probe conductances as a function of gate 
voltage G{V g ) for Fabry-Perot cavities on the single- and bi-layer devices, 
respectively, at zero magnetic field. We shift the curves horizontally by an 
amount V D which we identify as the gate voltage at which the Fermi level lies 
closest to the Dirac point. The conductance rises away from V g - V D = as 
observed by previously 2 ' 3 ' 4 . Small reproducible fluctuations of the conductance 
with magnitude of order e 2 lh can be seen (see insets to Figs. 2E and 2F); these 
fluctuations are not universal conductance fluctuations (UCF; see Supplemental 
Information) but, as argued below, result from the interference of ballistic electron 
waves in the Fabry-Perot cavity \ 

Figure 3 shows color-scale maps of the differential conductance d//d\/as 
a function of bias voltage V applied between the two electrical contacts and gate 
voltage V g for the single- and bi-layer devices shown in Figure 2. A pattern of 
diagonal lines of increased conductance is evident; this pattern is the signature of 
Fabry-Perot interference in a mesoscopic device 1,6 . Neighboring diagonal lines 



have similar slopes, and diagonal line of similar positive and negative slope are 
found in each V g region. Each individual diagonal line results from the 
enhancement in conductance when a particle-in-a-box resonance, or a group of 
constructively-interfering resonances, is aligned with the source electrode (+V) or 
drain electrode {-V); the symmetry about V = reflects the source-drain 
symmetry of the device. Note that the pattern is inconsistent with Coulomb 
blockade; there are no diamond-shaped low conductance regions around V = 0, 
and the overall conductance >> e 2 /h excludes Coulomb blockade. 

Resonant transmission through a Fabry-Perot cavity has been reported 
previously for carbon nanotubes (CNTs) 6,11 and graphene 1 . In the case of 
CNTs, there is a single path length L connecting the electrodes, and the 
resonances are evenly spaced in V and V g (see Methods). In graphene \ the 
resonances are randomly spaced, which may result from a spread of path 
lengths due to non-parallel electrodes or electron paths which are not 
perpendicular to the electrode-graphene interfaces. However, important 
information can be gained by analyzing the slope AV/AV g of the resonant lines in 
Fig. 3. Briefly, the slope measures the change in energy AE = eAV/2 (the factor 
of two results from the potentials +AV/2 and -AV/2 applied to the two electrodes 
relative to the graphene in a ballistic device) of the resonance as the particle 
number dn = c g AV g /e is changed. The slope is then equal to AV/AV g = 
(2c g /e 2 )AE/An = (2c g /e 2 )D" 1 ; i.e. the slope is inversely proportional to the density 
of states (see Methods for a more rigorous derivation of the same result). From 
Eqns. 3a and 3b above, we expect that AV/AV g ~ n' V2 ~ \V g - V D \' V2 for the 



single-layer (massless dispersion) sample, and AV/AV g ~ constant for the bi- 
layer (massive dispersion) sample. Figs. 3A-C show that the slope indeed varies 
significantly with gate voltage (electron density) for the single-layer graphene 
sample, with the highest slope occurring near V g - V D = 0. The slope is nearly 
constant in the bi-layer graphene sample (Figs. 3D-F). 

Fig. 4 plots the density of states D = (e 2 /2c g )(AV/A\/g)" 1 for the single- and 
bi-layer graphene samples extracted from Figs. 3A-F and additional data (not 
shown) as a function of electron density n = c g {V g - V D )/e. Data from an 
additional single-layer sample is also shown (see Supplementary Information). 
Solid lines are fits to Eqns. 3a and 3b for the single-layer and bi-layer data 
respectively. The expected dependences on particle number are verified: D ~ 
n V2 for massless particles in single-layer graphene (Eqn. 3a), and D ~ constant 
for massive particles in bi-layer graphene (Eqn. 3b). Only a single fitting 
parameter is used in each fit, v F =(1.09+0.01)xl0 6 m/5 for massless particles in 
single-layer graphene and m* = (0.0315±0.0001)m e for massive particles in bi- 
layer graphene. As discussed in detail below, the parameters are in excellent 
agreement with theoretical and other experimental results. 

We now discuss the detailed dependence of the density of states on 
particle number in single- and bi-layer graphene, and the implications of the 
results for understanding the electronic structure of these materials. From the fit 
to Eqn. 3a in Figure 4, we determine a Fermi velocity for single-layer graphene of 
v F =(1.09 + 0.01)xl0 6 m/s . A tight-binding model of graphene 12 gives 

Vp =(V3/2W /r| ~ 1.0 x 10 6 m/s where a = 2.46 A is the graphene lattice 



constant and Yo ~ 3.16 eV is the nearest-neighbor hopping parameter. The 
inclusion of electron-electron interactions will renormalize the Fermi velocity 
slightly 13 , and the slightly higher v F observed here is consistent with other 
experiments on graphene 2,3,14 ,15 . 

The density of states in single-layer graphene remains finite as n — > due 
to charge inhomogeneity caused by charged impurities near the graphene, as 
has been observed previously. The minimum density of states D on order 2 x 
10 12 eV" 1 cm" 2 corresponds to a charge density n on order 10 11 cm" 2 , in agreement 
with theoretical 16 and experimental expectations 8,17 for the minimum charge 
density at the Dirac point in the presence of charged impurity disorder due to the 
SiC>2 substrate. 

We now discuss bilayer graphene. From the fit to Eqn. 3b in Figure 4, we 
determine m = 0.032m e . Assuming v F = 1.09 x10 6 m/s we have yi = 0.40 eV, in 
excellent agreement with the experimental values for graphite of 0.39 ± 0.01 eV 
18 and with other experiments on bilayer graphene 19, 20 . Because the bands are 
not strictly parabolic, the density of states should depend on particle density, 
increasing with increasing particle density. The hyperbolic nature of the bands 
becomes important for particle densities roughly greater than y 1 2 /(4tt/7 2 v f 2 ) ~ 3 x 
10 12 cm" 2 . Experimentally, we see little variation in the density of states for 
particle densities up to 6 x 10 12 cm" 2 , indicating a wider range of validity of the 
parabolic spectrum than expected. We do not currently understand this 
discrepancy, but we note that electron-electron interactions should again be 
important, as was pointed out previously in the failure of the single-particle 



picture to quantitatively explain the cyclotron resonance spectrum in bilayer 
graphene 20 . 

In conclusion, we have probed the density of particle-in-a-box states as a 
function of particle number for massless 2d Fermions (single-layer graphene) 
and massive 2d Fermions (bi-layer graphene) in phase-coherent measurement. 
Understanding of coherent transport is an essential step to realize other 
interesting experiments in graphene such as focusing lens, klein tunneling and 
graphene superlattice 25,26 ' 27,28 . The density of states varies as the square-root of 
particle number for massless 2d Fermions, and is constant for massive 2d 
Fermions. The single parameters in the dispersion relations are extracted; the 
Fermi velocity v F =(1.09+0.01)xl0 6 m/s for massless particles in single-layer 
graphene and m* = (0.03 15 ± 0.000 \)m e for massive particles in bi-layer graphene, 

in excellent agreement with theoretical expectations and other experimental 
observations. 
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Figure Captions 

Fig. 1. (A) Massless dispersion and (B) massive dispersion relations in two 
dimensions. (C-D) Allowed wavenumbers for particle in a box of aspect ratio W/L 
= 1 .6. Solid lines are contours of equal energy for massless dispersion relation 
(C) and massive dispersion relation (D). (E-F) Particle energy as a function of 
particle number in a box with W/L = 1 .6 for massless dispersion relation (E) and 
for massive dispersion relation (F). 

Fig. 2. (A-B) Optical micrographs of single-layer graphene device (A) and bilayer 
graphene device (B). (C-D) Longitudinal and Hall conductivity as a function of 
gate voltage at magnetic field of 9 T and temperature of 1 .3 K for single-layer 
graphene device (C) and bilayer graphene device (D). (E-F) Two-probe 
conductance as a function of gate voltage at zero magnetic field and temperature 
of 1 .3 K for single-layer graphene device (E) and bilayer graphene device (E) 

Fig. 3. Color-scale two-dimensional plots of differential conductance G = d//dVas 
a function of bias voltage Vand gate voltage V g measured in single (A-C) and 
bilayer (D-F) graphene at temperature T = 1.3K. A smooth background 
conductance was subtracted to enhance the patterns. The sample dimensions 
are 1 .5um(l/V) x 0.3um(/_) for single layer graphene and 4.3um(l/V) x 0.2um(/_) for 
bilayer graphene. Yellow lines illustrate the slope of the Fabry-Perot resonances. 



Fig. 4. Density of states of single-layer graphene (blue symbols; data from two 
devices shown) and bilayer graphene (red symbols, data from one device shown) 
as a function of particle density. Solid lines are fits to Eqn. 3a (blue) with v? = 
1 .09 x 1 6 m/s and Eqn. 3b (red) with m* = 0.032 m e . 
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